astro-ph/0106022 



Oscillations During Inflation 
and 

the Cosmological Density Perturbations 



N. Bartolo^ 2 S. Matarrese^ ^ g^^^^ Riotto^ 

^ Dipartimento di Fisica di Padova "G. Galilei" 
Via Marzolo 8, Padova 1-35131, Italy 

and 

"^INFN, Sezione di Padova 
Via Marzolo 8, Padova 1-35131, Italy 



Abstract 

Adiabatic (curvature) perturbations are produced during a period of cos- 
mological inflation that is driven by a single scalar field, the inflaton. On 
particle physics grounds - though - it is natural to expect that this scalar 
field is coupled to other scalar degrees of freedom. This gives rise to oscil- 
lations between the perturbation of the inflaton field and the perturbations 
of the other scalar degrees of freedom, similar to the phenomenon of neu- 
trino oscillations. Since the degree of the mixing is governed by the squared 
mass matrix of the scalar fields, the oscillations can occur even if the energy 
density of the extra scalar fields is much smaller than the energy density of 
the inflaton field. The probability of oscillation is resonantly amplified when 
perturbations cross the horizon and the perturbations in the inflaton field 
may disappear at horizon crossing giving rise to perturbations in scalar fields 
other than the inflaton. Adiabatic and isocurvature perturbations are in- 
evitably correlated at the end of inflation and we provide a simple expression 
for the cross-correlation in terms of the slow-roll parameters. 
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1 Introduction 



It is commonly believed that the Universe underwent an early era of cosmological in- 
flation. The flatness and the horizon problems of the standard big bang cosmology are 
elegantly solved if, during the evolution of the early Universe, the energy density is 
dominated by some vacuum energy and comoving scales grow quasi-exponentially. The 
prediction of the simplest models of inflation is a flat Universe, i.e. Qtot = 1 with great 
precision. 

Inflation [|l| has also become the dominant paradigm for understanding the ini- 
tial conditions for structure formation and for Cosmic Microwave Background (CMB) 
anisotropy generation. In the inflationary picture, primordial density and gravity-wave 
fluctuations are created from quantum fluctuations "redshifted" out of the horizon, 
where they are "frozen" as perturbations in the background metric ^ ^, |, Q . Metric 
perturbations at the surface of last scattering are observable as temperature anisotropy 
in the CMB. 

Primordial perturbations can be of two kinds, the adiabatic and the isocurvature ones. 
Recently a lot of attention has been drawn on correlated mixtures of the two 0, |^. 
In particular in Ref. [§] it has been shown how the correlation between the adiabatic 
and the isocurvature mode gives rise to new features both in the CMB anisotropics and 
in the large scale structure. These scenarios are strongly different from those usually 
consid ered up to now, in which only independent mixtures of the two modes were con- 
sidered. 

Analyses of the possible constraints on these perturbations coming from the present 



data 1^ and future experiments [10| have been also made, as well as investigations on 
the production of the correlation during inflatio n: in Ref. a speciflc model of double 



inflation was considered, and in Ref. |TT] a transparent formalism for studying the adia- 
batic and isocurvature modes was introduced. Here we pursue further the investigation 
on the possible correlation mechanism, trying to be as general as possible. 
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Our starting point is the simple observation that, if the inflaton field couples to other 
scalar degrees of freedom, oscillations between the perturbation of the infiaton field and 
the perturbations of the other scalar degrees of freedom are induced. For this phe- 
nomenon to happen, it is sufficient that the mass squared matrix of the scalar degrees 
of freedom is not diagonal. This induces a mixing among the different scalar states and 
such a mixing can be large even if the energy density of the inflaton fleld dominates over 
the energy density of the other scalars. We will show that the probability of oscillation is 
resonantly amplifled when perturbations cross the horizon. Adiabatic and isocurvature 
perturbations are therefore correlated at the end of inflation and we provide a simple 
expression for the cross-correlation in terms of the slow-roll parameters. 

The plan of the paper is as follows. In Section 2 we show that, during an inflationary 
period in which several scalar flelds are present, is natural to expect a mixing and 
consequent oscillations between the fluctuations of the scalar flelds. In Sectio n 3 the 
correlation of the adiabatic and isocurvature perturbations is explained in terms of this 
oscillation mechanism, and an explicit expression for it is derived. Finally, in Section 4, 
we show how to set the initial conditions for structure formation in the post-inflationary 
epoch in the case where the isocurvature perturbations and the correlation are present. 

2 Oscillations during inflation 

Despite the simplicity of the inflationary paradigm, the number of inflation models that 



have been proposed in the literature is enormous |]12| . Models have been invented which 



predict non-Gaussian density fluctuations |jT3[, isocurvature fluctuation modes fT^, and 
cosmic strings |]T5| . 



The simplest possibility is represented by the so-called single-fleld models of infla- 
tion, where the the inflationary epoch can be described by a single dynamical order 
parameter, the inflaton fleld 0. Quantum fluctuations of the inflaton fleld produce 
Gaussian adiabatic perturbations of the metric with a nearly scale independent spec- 
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trum, ns — 1. The amplitude of the perturbation can be characterized by the co moving 
curvature perturbation TZ, which remains constant on super-Hubble scales until the 
perturbation comes back within the Hubble scale long after inflation has ended. Single 
models of inflation have already been started to be constrained by the recent accurate 
measurements of the CMB anisotropy |17| |. 

On particle physics grounds - though - it is hard to believe that only one single 
scalar field plays a role during the inflationary stage. On the contrary, it is quite 
natural to expect that during the inflationary dynamics several other scalar fields xi 
(J = 1, ■ ■ ■ , A^) are present. As soon as one considers more than one scalar field, one 



must also consider the role of isocurvature fluctuations [T8[ . Such perturbations produce 
an anisotropy in the Cosmic Microwave Background radiation which is six times larger 
than the adiabatic perturbations. To obtain small CMB anisotropy and still explain 
galaxy formation one has to strongly suppress isocurvature perturbations on the horizon 
scale while keeping them sufficiently large on galactic scales. This can be done if the 
additional scalar fields acquire a mass of the order of the Hubble rate H during inflation. 
Furthermore, the presence of more than one scalar field may not only affect the evolution 
of the curvature perturbation, but also give rise to the possibility of seeding isocurvature 
perturbations after inflation. 

The contribution to the total energy density of the extra scalar degrees of freedom 
Xi might or might not be negligible compared to the one provided by the scalar field 
0. If it is, the model of inflation is called multiple field model, a general formalism to 
evaluate the curvature perturbation at the end of inflation in such models was developed 



in Ref. 19 



However, even if the contribution to the total energy density of the scalar fields xi 
is small, quantum fluctuations 6xi of the scalar fields are amplified by gravitational 
effects. To compute the amplitude and the spectral shape of such perturbations one 
may consider the theory of a single free scalar field 6xi with mass m in a de Sitter 
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background . 

In this section we wish to show that the generation of the quantum fluctuations in 
the fields xi may be due to another mechanism, which we call oscillation mechanism 
and it is important to stress that such a novel mechanism operates even if the energy 
density of the scalar field xi is much smaller than the contribution coming from a single 
scalar field (p which - for such a reason - deserves the name of infiaton field. 

If the scalar sector is composed by more than one single scalar field 0, all the scalar 
degrees of freedom will in general mix. In the particle physics language this can be 
translated by saying that the mass squared (or the Hamiltonian) in the basis (0, xi) 
is not diagonal or, equivalently, that the states (0, xi) are interaction eigenstates, but 
not mass eigenstates. If the Hamiltonian of two quantum states is not diagonal, the 
interactions eigenstates (0, xi) oscillate during the time evolution of the system. A 
familiar example in particle physics is represented by the oscillating system of kaons 
and K°. 

The oscillation mechanism responsible for the amplification of the quantum fluctu- 
ations of the fields xi is due to the oscillations of the infiaton fluctuations into fluctua- 
tions of the scalar fields xi^ iii the presence of the infiaton background field 0o. In other 
words, a perturbation in the infiaton field may evolve (oscillate) into a perturbation 
of another scalar degree of freedom x with a calculable probability. 

From a pure quantum mechanical point of view, such fluctuations are generated 
as coherent states. Therefore the oscillation mechanism is a coherent production of 
quantum fluctuation^. These oscillations during inflation show a behaviour similar to 
the one present in the phenomenon of coherent neutrino oscillations in a medium like 
the Sun or the Earth. As we shall see, oscillations during inflation are charecterized by 
amplification effects analogous to the MSW effect of solar neutrinos. 

Before launching ourselves into the details, let us now give a couple of examples 



-'^The importance of coherent production of particles after inflation during the preheating stage has 
been recently emphasized in Ref. pl[. 
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supporting the fact that a large mixing between different scalar fields is generically 
expected during the inflationary stage. In supergravity and (super) string models there 
exists a plethora of scalar fields - loosely called moduli - with gravitational-strength 
couplings to ordinary matter. The mass of these scalar fields is of the order of the 
gravitino mass m^j2 ~ 100 GeV in the present vacuum, but is of the order of the Hubble 
rate H during inflation. In these theories coupling constants and masses appearing in 
the Lagrangian have to be thought as functions of the dimensionless ratio x//^Pb 
where xi (-^ = ' ' ' 5 ^) denotes a generic modulus field and Mpi is the Planck mass. 
For instance, a generic coupling constant A is in fact a function of the scalar moduli 

where c is a coefficient usually of order unity and 5xi = Xi ~ ixi)- This expansion 
introduces a direct coupling between the scalar moduli and the infiaton. The potential 
V becomes a function of two (or more) fields, V = V{(j), xi) and the second derivative 
{d'^V/d(j)dxi) may be as large as during inflation. If we set = Xoy all elements 
of the mass squared matrix Aifj = {d'^V /dxid(t)j) {j = 0, ■ ■ ■ , N) are of the form Cij H"^ 
where Cij = 0{1). A considerable mixing between the infiaton field and the moduli 
fields is generated if the infiaton background field 0o takes values as large as the Planck 
scale. Notice that this may happen even if cqo <^ 1, as required by the flatness of the 
infiaton potential. Under these circumstances, the perturbation of the scalar field 
may oscillate into a perturbation of the modulus field which is generated as a coherent 
state. 

Another example may be provided by theories in which gravity may propagate in 
extra-dimensions p2| where there appears a infinite tower of spin-0 graviscalar Kaluza- 



Klein excitations and the infiaton field can mix to these particles by coupling to the 
higher-dimensional Ricci scalar. 

Let us now describe a simple, but illustrative example. Consider two scalar fields, 
and X- We will dub cf) the infiaton field, even if this might be a misnomer as the two 



fields might give a comparable contribution to the total energy density of the Universe. 

The scalar field perturbations, with comoving wavenumber k = 27ra/A for a mode 
with physical wavelength A, obey the perturbation equations 



^2 

5(f) + 3H5(l) + —S(f) + V^^5(l) + V^Jx = 
a 

Sx + 3H5x + ^Sx + V^Jx + V^4,5cl> = 0, (2) 

where we have indicated by V^,/, = {d^V / d(j)d(j)) and similar notation for the other 
derivatives. 

The squared mass matrix is given by 

We now introduce a time-dependent 2x2 unitary matrix U such that 

U'^M'^U = diag {ujI, ujI) = uj'^. (4) 

In the following we will assume that all the entries of the squared mass matrix M.'^ are 
real, so that the unitary matrix U reduces to an orthogonal matrix 

^'^^^ — sin^ \ , s 

"~ sin^ cos^ ' 



where 



tan2e= /\ (6) 



and the mass eigenvalues are given by 

1 



a;?n = - 

1,2 2 



(^^<A + ^xx)±V(n0-^xx) +4K 



x4> 



(7) 



Adopting the vectorial notation ^ — (^i, ^2)^ = x)^, the equations of the scalar 

perturbations may be rewritten as 



S-^ + (3H + 2U^U) 5^ + + a;2 + 3HU'^U + U^U^ 5^ = 0, 



(8) 



where 

ijTfj = ( * 

-1 

and 



U^u^i \ \ ]e (9) 



u-u^[\\)u['^ \ )e^. (10) 

Notice that the matrix U diagonahzes the squared mass matrix at the price of intro- 
ducing further non-diagonal terms in the equation of motion. 

To proceed further, we may now take advantage of the slow-roll conditions which 
are to be attained during inflation lO]. If we consider the generic slow-roll parameters 



= r,., = M^^, (11) 

where M = Mpi/-\/87r is the reduced Planck mass, a successfuU period of inflation 
requires that \eij,riij\ <C 1, i.e. the potential has to be flat enough for inflation to 
develop. 

Since time derivatives of the slow-roll parameters are second-order in the slow-roll 
parameters themselves, ?) ~ (^{e^, rf ), it is easy to convince oneself that - if we keep 
only the slow-roll paramaters at the first-order - Eq. (^ gets simplified to 

5^ + ?,H5i! + [^+uj'^ + 3HU'^l?j5^ = Q. (12) 

Introducing the conformal time dr = dt/a, where a is the scale factor of the expanding 
Universe, and the rescaled fields ^'^'i = aS'^i and 6^2 = Eq- (HID becomes 

^ a" ^ 



S^'l +{^^-^+ ] + 3^ O'a 6^2 = 



6^2+ (k'^ -^+^Wj^'^2-^He'a6^i = (13) 

To illustrate the phenomenon of oscillations during inflation, we make the assump- 
tion that the non-diagonal terms in Eq. (|13D proportional to 3H6'a = 3H6a^ are 
smaller than the diagonal entries cj^g'^^- This hypothesis is correct, for instance, if the 
squared masses are constant in time. We adopt this simplification in order to render 
the description of the phenomenon of oscillations during inflation more transparent. 



In fact, in the following section we will solve the problem exactly at the first-order in 
the slow-roll parameters and show that at this order of approximation taking 3H6'a <^ 
uf 2 does not change the final result for the probability. 

Eq. (|T3]) is solved byQ 

S^, = l0Fe<'^"+^)5 {-rY/'Hj^\-kT), z = 1,2, (14) 

where the conformal time r assumes negative values (the beginning of inflation is at 
some |r| ^ 1), Hj^^^ are the Hankel's functions of the first kind and uf = 9/4 — (cjj/if)^. 



The normalization factor in Eq. (|l^) is chosen such that matches the plane-wave 
e~'^'''^ / V2k for subhorizon scales in the far ultraviolet k/aH ^ 1. 

We are now in the position to ask what is the probability (as a function of time) 
that a scalar perturbation in the "infiaton" field becomes a scalar field perturbation 
in the scalar field 6x- The answer to this question is readily given if we remember 
that [0(r),x(r)]^ = W[^i(r), ^2(t)]^- This means that - at a given time r - the 
scalar perturbations 6(j) and 6x are a linear combination of the mass eigenstates scalar 
perturbations 5\E'i and 6"^ 2 

50=5:^1,5^,, 6x=J2^2iS^e. (15) 

i=l,2 1=1,2 

The probability that a scalar perturbation 5(f) at the time tq becomes a scalar pertur- 
bation 5x at the time r is therefore given in general by 



P[5<P{t,)^6x{t)] 



?=i,2 \ \o^i\j) 



(16) 



We take as initial condition tq — 00 , and we follow the conversion probability at 
time r. 

If the unitary matrix U is real, Eq. (ITB]) becomes 



^To be consistent one should also expand the factor ^ up to the first order in the slow-roll pa- 
rameters, reflecting the fact that inflation does not generically occurs with a pure de Sitter dynamics. 
However, these corrections do not change the final expression for the oscillation probability, see Eq. 
(pO|), since they would alter the wave- functions of the two mass eigenstates in an equal manner. 
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Let us investigate how such probabihty changes as a function of the physical wavelength 
A = 2Ti{a/k). At subhorizon scales, k > aH, the functions tend to the common 
plane-wave solution e~^^'^ / \/2k and therefore 

P [50(ro) ^ 5x{r)] ~ {k->aH). (18) 

However, at superhorizon scales, k < aH, the functions 5'^i develops a phase dependence 

6^1 1 ~ e*('^^-5)f 2"^-i ^ (-A;r)5-'^* (19) 

r(3/2) ^ ^ ^ ^ 

and the conversion probability becomes 



1 — cos ( — Az/ 
2 



sin^ 2^ sin^ ( ^ Az/ ) {k<^aH), (20) 



where Az/ = z/i — z/2. 

In the limit cu^g such a probability reduces to 



TT Acj2^ 



P [50(ro) ^ 5x(r)] =^ sin^ 29 sin^ ( ) , (21) 



with Acj^ = ujI — uj. 



Note that the same formula holds for the probability P [(5x(''o) ^'Pi'^)]- This is 
due to the fact that eq. ([16|) depends only on the mass eigenstates and it is symmetric 
in 5§i and 5§2- 

The expression (|20|) reminds the well-known formula which describes the evolution 
in time of the probability of oscillations between two neutrino flavours p3|. In both 
cases the probability is identically zero if the two mass eigenstates have equal masses 
(no oscillations are present in the degenerate case); there is the same dependence (as 
sin^ 26) on the mixing angle 6 and the same functional dependence on the difference of 
the squared masses. Differences are present, though. While the probability of neutrino 
conversion is depending upon time for any value of the neutrino energy, in our case at 
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superhorizon scales the conversion probability becomes constant in time. This does not 
come as a surprise since on superhorizon scales the dynamics of the system is frozen. 

What is more interesting is the time evolution of the conversion probability. Consider 
a scalar perturbation in the inflaton field with a given physical wavelength A = 2TT{a/k) 
which gets stretched during inflation. As long as the wavelength remains subhorizon, 
the scalar perturbation remains a pure perturbation in the inflaton field. However, 
as soon as the wavelength crosses the horizon, the perturbation in the inflaton field 
may become (generate) a perturbation in the other scalar field x with a nonvanishing 
probability determined by Eq. (|20|). At horizon crossing there is an amplification 
mechanism of the fluctuations in the field x which is reminiscent of the MSW effect 
operative for solar neutrinos. 

The phenomenon of resonant amplification is easily understood if one remembers 
that a given wavelength crosses the horizon when k = aH, i.e. when k"^ = a" /a using 
the conformal time. As long as the wavelength is subhorizon, k"^ ^ a"/a, the presence 
of the mass terms in the equations of motion (^) is completely negligible compared to 
the factor {k"^ — a" /a). On the other hand, when the wavelength crosses the horizon 
the term [k"^ — a" /a) vanishes and the effect of the mixing in the mass squared matrix 
is magnified, giving rise to the resonant effect. Finally, when the wavelength is larger 
than the horizon, k"^ <^ a" /a, the term {k"^ — a" /a) starts to dominate again over the 
mass terms and the oscillations get frozen. 

We conclude that fluctuations in the scalar field x ^"^^ generated as coherent states 
at horizon crossing through the oscillation mechanism out of perturbations in the scalar 
field (j). 

A couple of comments are in order here. First of all, we wish to stress that the 
oscillation mechanism operates even if the energy of the inflaton field is much larger 
than the energy stored in the other scalar field x- This is because what is crucial for 
the oscillations to occur is the relative magnitude of the elements of the mass squared 



10 



matrice M?' . Secondly, the magnitude of the probabihty depends upon two quantities, 
sin^ 1Q and Acu^/ iP'. Both can be readily expressed in terms of the slow-roll parameters. 
The first factor is not necessarily small, in fact it may be even of order unity for maximal 
mixing. If expanded in terms of the slow-roll parameters, it is 0{rf ^ e°) . The second term 
is naturally smaller than unity and is linear in the slow-roll parameters. This reflects the 
fact that during inflation only perturbations in those scalar fields with masses smaller 
than the Hubble rate may be excited. However, ISuP- jH"^ is not necessarily much smaller 
than unity and the amplification of the conversion probability at horizon crossing may 
be sizeable. 

3 Correlation between adiabatic and entropy per- 
turbations during inflation 

As we already mentioned in the previous section, adiabatic (curvature) and entropy 
(iso curvature) perturbations are produced during a period of cosmological infiation if 
more than one scalar field is present at this epoch. In this section we wish to provide a 
simple expression for the cross-correlation between adiabatic and entropy perturbations 
inspired by the considerations developed in the previous section. 

A consistent study of the hnear field fiuctuations requires the knowledge of the linear 
scalar perturbations of the metric, corresponding to the line element 

+ a^{{\-2i))bij^2E^^\dx'dxK (22) 

The consistent equation for a generic scalar field perturbation ^xi = ' ' i^) 
with comoving wavenumber k — 27ra/A for a mode with physical wavelength A reads 

^Xi + ^H5xi + —^Xi + ^xixJXj 



0? J 



= -2V^,A + X: 



A + 3ip + —(a'^E-aB) 



(23) 
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At this stage it is useful to introduce the gauge-invariant Sasaki- Mukhanov variables 

m 



Qi = 5ipi + (24) 
which, in the spatially flat gauge = 0, obey the equations of motion 



Qi + 'iHQi + —Qi + M]jQj = 0, (25) 



where 



^iJ = ^x.x. - ^ hXiXj - ^ (^r,.. - 3 e.. ) , (26) 



where the last expression has been obtained performing an expansion up to the first 
order in the slow-roll parameters. This equation is similar to Eq. (Q) and from what we 
have learned in the previous section, oscillations among the different quantities Qj are 
expected to take place. Notice also that the Q/'s oscillate even if the part of the mass 
squared matrix A^/j proportional to V^x/x.j diagonal. This is because non-diagonal 
entries are always present because of the e/j-parameters. 

Once the variables Qj have been defined, one can define the comoving curvature 
perturbation p5|, ^ 



and give a dimensionless definition of the total entropy perturbation (automatically 
gauge- invariant ) 

5 = i/f^-^V (28) 
\P P J 

For scalar fields the latter is given by 

Let us now restrict ourselves to the case of two fields, and x- Following the nice 



treatement of Ref. [11|, we can define two new adiabatic and entropy fields by a rotation 



in field space. We define the "entropy field" s |11 



Ss = (cos (3)Sx - (sin (3)6^, (30) 
12 



where 



cos [3 



I V,2 



sin j3 



Ji 

\2 ^2 



(31) 



Notice that bs can be rewritten as 



= (cos/5)(5x - {sin /3)Q^, 



(32) 



From this definition it follows that s =constant along the classical trajectory, and hence 



entropy perturbations are automatically gauge-invariant [p7| [Till . 

The adiabatic part of the perturbation is associated to the orthogonal combination 



5Qa = {sin i3)Q^ + {cos i3)Q^. 



(33) 



Our goal is to give an expression of the cross- correlation between the adiabatic and the 
entropy perturbations 



{Q^{k)5s*{k')) = ^CQ^Ss5{k - k') 



(34) 



To do so, we adopt the technique developed in the previous section. As we have seen, 
though, introducing a unitary matrix U which diagonalizes the mass squared matrix 
is not enough to diagonalize the full system, (see Eq. (p^). This happens because, 
in general, the fields are coupled together. To proceed, we first define the comoving 
fields Qfj) = aQtf) and = aQ^, then we introduce a basis for annihilation and creation 
operators and a| (i = 1, 2) and perform the decomposition (r is the conformal time) 



Qx 



Q<t, 




ikx 



ikx 



h{r) 



h{r) 



ai{k) 
a2{k) 

ai{k) 
a2{k) 



+ h.c. 
+ h.c. 



(35) 



where 11^^ and Ilg are the conjugate momenta of Q^j, and respectively, and h and 
h are two 2x2 matrices satisfying the relation 



h h* - h* 
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i 6, 



(36) 



derived from the canonical quantization condition. The matrix U is given in Eq. (^. 
As it diagonahzes the squared mass matrix ( p6D (with the identification Xi = and 
X2 = X)i the mixing angle is given by 



tan 26* 

and the mass eigenvalues are given by 

1 



(37) 



(3^ 



It is not difficult to see that the matrix h satisfies the following differential equation 



h" + 2UU' h' 



/i = 0, 



(39) 



If we now expand up to the first order of perturbation in the slow-roll parameters, we 
obtain 



/ill + 2 ^' /i2i + + ]hu + H 9' a /lai 



a 



/I'a'i - 2 0' /I'li + ( - ^ + ulc? ] - H 9' a 



0, 



0, 



(40) 



and similar equations for hi2 and with the substitution hn hi2 and /i2i /i22- 
Making use of the decomposition (|35|) and definitions (|3^ ) and (p3D, we find that 



a^{QA{k)5s*{k')) = {sfsce - Cfsse) {c(}Ce + s,3Sg) 

+ {cpce + spsef [hnhl^ + h^h 



[spce - cpsef [/i2i/in + h22hl2] 



\h22?' - \hii? + \h2i? - \hi2?' 



(41) 



where we have made use of the shorthand notation C/3(e) = cos (3{9) and s^i^e) = sin P{9). 
To check this expression, we can consider the following limiting cases: 

i) No squared mass matrix: If the squared mass matrix (p6|) is identically zero, 
and are already orthogonal states and mass eigenstates and their time evolution 
is identical. This means that there is no correlation between them and one expects 
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{QA{k)5s*{k')) oc spcp ({Q^{k)Ql{k')) - {Q^{k)Q*^{k'))) = 0. This is confirmed by Eq. 
(PI) since in this case cg = 1, sg = 0, hu = /i22 and hi2 = /i2i = 0. 

ii) Two noninteracting fields: Suppose the system is composed by two fields (p and 
X with potential V = im^(0^ + x^)- In such a case, cp = Sjs = 1/V2 and the squared 
mass matrix (|2^) as the following structure 

-M^«(} J), (42) 

coming from the terms proportional to the e-parameters. This means that the mixing 
angle 9 = 7r/4 and cg = sg = 1/v^. The system is already diagonalized by the 
unitary matrix U and there is no need for the nondiagonal elements of the matrix h, 
hi2 = = 0. In such a case, the expression (|l|) predicts that {QA{k)5s*{k')) vanishes. 
This confirms the findings of Ref . . 

To compute the cross-correlation between the adiabatic and the isocurvature modes 
in a more general way, we can expand the solutions of the system (|40|) in powers of the 
slow-roll parameters. 

In the far ultraviolet k ^ aH the squared mass matrix is negligible and and 
are mass eigenstates. Therefore, we set the physical initial conditions by posing the 
unitary matrix lA proportional to the unity matrix and hij = 5ij e~^^'^ / \/2k. A simple 
inspection of the system (^) tells us that the functions hi2 and h2i are sourced by 
the functions /i22 and —hu, respectively. Since 6' is already 0{rj,e) during the time 
evolution of the system hi2 = — /i2i = C^iv^^)- Expanding Eq. (|4l|) up to first order in 
the slow-roll parameters, we get 

a^{QA{k)Ss*{k')) = {spCg - CpSg) {cpCg + SfsSg) 

+ [hiih*2, - h2ihl,] , (43) 
where we have made use of the fact that \h22f — \huf = 0{ri, e). 



h22\ - I /ill 
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A simple perturbative procedure shows that 

/(or(r)-r(o/(r) 



h2i(r) 



(44) 



where 



fir) = ^-V^e^-i-ry/'Hi%{-kr), 



w{t) = fr-rf 



(45) 



h 



22 



-ikr 



/V2k 



As long as the wavelength of the perturbations is subhorizon, hu 
and solving Eq. (|4^) gives |/i2i| oc (H/k) \hii\ <^ 

At superhorizon scales k <^ aH, hu, /i22 are proportional to the scale factor (up to 
terms linear in the slow- roll parameters). Since 9' = aO = a H F{e,T]), where F is a 
function linear in the slow-roll parameters, the integral in Eq. (0) can be performed 
and gives 

/^2i = -3/iiiF(e,r/)ln(^-^^. (46) 

With this information at hand, let us first evaluate the probability that a the per- 
turbation Qif, at a time tq becomes a perturbation at a generic time r. Using the 
decomposition (|35|) and making use of the fact that hi2 = — /i2i = 0{ri, e) during the 
time evolution of the system, the probability may be written as 



P[Qdro)^QAr)] = \{Qxir)mro))\' 
hii{ro) K-^ir) /i22(ro) /i22(^) 



ce se 



\hu{ro) h*u{T)\ \h22{To) h*22{T)\ _ 

Such a probability is vanishing at subhorizon scales, but at superhorizon scale 

1 



(47) 



sin^ 29 



1 — cos ( —Au 
2 



sin^ 29 sin^ [-rAu] , 



TT 



(4J 



which reproduces Eq. (pO|). 

Therefore, we expect that the cross-correlation is extremely small at subhorizon 
scales reflecting the fact that and are good mass eigenstates Nevertheless, the 
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cross-correlation does not vanish at superhorizon scales and is given by 



Cqj^ss = {spce - cpse) {cpce + spse) ( ^ j 



1 - 



~aH 



(49) 



where Hk — k/a and Av — V2—V1 (which is given by — 1/3 Aa;^/if^ = —1/2){ujI—ujI)/H'^ 
if a;^ 2 <^ H^)- If we normahze the scale factor a in such a way that a = 1 at the end of 
inflation and indicate by iV = —In a the number of e-folds until the end of inflation, we 
finally get 



1- - e3 ^ 



(50) 



This result is remarkably simple and can be easily expressed in terms of the slow-roll 
parameters using the expressions 



~ 3 



''l(p(t> Vxx 2 i^4>4> ^xx) 



+ 4 {r}4>x - ^Hx 



and 



tan 2^ = 



2 (^</.x - \Hx) 



Vxx 3 (^'t"t> ^xx) 



(51) 



(52) 



The origin of the cross-correlation is due to a rather transparent physical behaviour. 
At the inflationary epoch, the gauge invariant perturbations and are generated 
with different wavelengths stretched by the superluminal expansion of the scale factor. 
Since the squared mass matrix of Q,^ and is not diagonal, oscillations between the 
two quantities are expected. Till the wavelength remains subhorizon, and Q^. evolve 
independently and may be considered good mass eigenstates. However, as soon as the 
wavelength crosses the horizon, an amplification in the probability of oscillation be- 
tween and occurs: a nonvanishing correlation between and is created on 
superhorizon scales because of the nondiagonal mass matrix Mjj. Since the adiabatic 
and the isocorvature modes are a linear combination of and Q^, at horizon crossing 
a nonvanishing correlation between the adiabatic and the isocurvature modes is left 
imprinted in the spectrum. 
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4 Initial conditions in the radiation era 



It must be emphasized that the mechanism described above for the production of adia- 
batic plus isocuvarture perturbations and their cross-correlation is active during infla- 
tion. A common feature of isocurvature perturbations is that they might not 

survive after the end of inflation [pSj |29| , |30| . If during reheating all the scalar fields 
decay into the same species (photons, neutrinos, cold dark matter and baryons), the 
only remaining perturbations will be of adiabatic t 

ype. 

In fact, in order to have isocurvature perturbations deep in the radiation era it is 
necessary to have at least one non-zero isocurvature perturbation |]3T| : 



Sc, ^ - ^ 0, (53) 

where 6a = Spa/pa , Wa = Val Pa (the ratio of the pressure to the energy density), and 
a and (3 stand for any two components of the system (they can also be scalar fields). 
Saji (a gauge-invariant quantity) measures the relative fluctuations between the different 
components. Adiabatic perturbations are characterized by having S^p = for all of the 
components. 

So in the following we will assume that the mixing between the scalar fields is negligible 
after inflation and that, for example, the scalar field decays into "ordinary" matter 
(photons, neutrinos and baryons) and the scalar field x decays only into cold dark 
matter (or x does not decay, like the axion, so that it consitutes the cold dark matter). 
In this case we can write: 

3 3 

^CDM = ScDM-rest + , = ^^-y = -^^u = h (54) 

where 5a specifies the amplitude of the adiabatic mode of perturbations, and "rest" 
stands fo photons, neutrinos and baryons. 

In fact the initial conditions for the evolution of cosmological perturbations are set 
once Sap and the gravitational potential $ (in the longitudinal gauge) are given deep 
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in the radiation era 0]. These initial conditions are necessary to explore the effects 
on the large scale structure of the universe and the CMB anisotropics from correlated 
adiabatic/isocurvature perturbations ||. 

Here we want to show how to link the value of Scdm -rest and $ to the inflationary 
quantities 6s and Qa- 

For the gravitational potential the main equation is: 

^ = -f + (l - ^) <^>, (55) 

where TZ is the curvature perturbation (see [O, and reference therein). 
Now, making an expansion in the slow-roll parameters to lowest order, it can be shown 
that the term proportional to $ can be neglected for most of the inflationary period. 
During the radiation dominated epoch —H'^/H = 1/2, and by matchi 
ng the two stages we can write: 

7^rad = (56) 

where 7?-rad is the curvature perturbation at the end of inflation, and it is directly related 



to Qa through ITT 



n = QA^ , i = (cos/3)0+ (sin/5)x. (57) 
As far as Sc dm -mst is concerned, let us define the following quantity: 

^x4> = — -^- (58) 
X (p 

For the two scalar fields (p and x the isocurvature perturbation S^^ as defined in Eq. 
(53) results to be S^^jj = a? {5^^/ a?) [^ . 
On the other hand: 

5s = ^ S^^ . (59) 

Then, to lowest order in the slow-roll parameters and taking into account only the 
growing isocurvature mode, one finds: 

S^^ = ^^^^ + 5s . (60) 
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To match to the radiation epoch we take ScoM-rcst = S^^j, at the end of inflation. 
Finally, we can give an expression for the correlation between $ and ScDM-rcst which 
represents the correlation between the adiabatic and isocurvature modes: 

{mS*{k')) = - ^^ {QA{k)5s*ik')), (61) 

where the right-hand side of this equation is evaluated at the end of inflation, and 



In a separate paper |3§] we calculate in detail the amplitudes and spectral indeces in 
terms of the slow-roll parameters for Sajs and the correlation. 



5 Conclusions 

In this paper we have shown that the correlation between adiabatic and entropy per- 



turbations during inflation 0, |T3] can be explained through an oscillation mechanism 
between the perturbations of two or more scalar fields. 

If the isocurvature perturbation mode survives after inflation, this correlation can leave 
peculiar imprints on the CMB anisotropics as an additional parameter for the initial 
conditions of structure formation. In the near future very accurate data on CM B 



anisotropics will be available from the MAP and Planck |^ experiments. Thus it 
is worth investigating further these issues, since they can represent a valid alternative 
to the simplest inflationary models. The latter are based on a single slow-rolling scalar 
field and generally predict adiabatic, scale-free and nearly Gaussian perturbations. 
In this respect the relevant point of this work is to stress how isocurvature perturba- 
tions and cross correlations come out in a very natural way even within a single-field 
inflationary scenario, with an inflaton (p which gives the dominant con tribution to the 
total energy density and other scalar fields whose energy densities may or may not be 
subdominant. Moreover, the same correlation can give rise to a transfer of possible 
non-Gaussian features from the isocurvature mode to the adiabatic one. Different infla- 
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tionary models, in fact, have been proposed in the hterature which predict non-Gaussian 
isocurvature perturbations ||13|. In the present case, however, even if the isocurvature 



component is suppressed, as the observations suggest 0, this transfer mechanism can be 
very efficient for producing non-Gaussian adiahatic perturbations. We will investigate 
this possibility in a separate paper 
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